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ABSTRACT 



Context. Optical high-resolution spectra of five central stars of planetary nebulae (CSPN) in the Galactic bulge have been obtained with 
Keck/HIRES in order to derive their parameters. Since the distance of the objects is quite well known, such a method has the advantage that 
stellar luminosities and masses can in principle be determined without relying on theoretical relations between both quantities. 
Aims. By alternatively combining the results of our spectroscopic investigation with evolutionary tracks, we obtain so-called spectroscopic 
distances, which can be compared with the known (average) distance of the bulge-CSPN. This offers the possibility to test the validity of model 
atmospheres and present date post-AGB evolution. 

Methods. We analyze optical H/He profiles of five Galactic bulge CSPN (plus one comparison object) by means of profile fitting based on state 
of the art non-LTE modeling tools, to constrain their basic atmospheric parameters (7^ ff , logg, helium abundance and wind strength). Masses 
and other stellar radius dependent quantities are obtained from both the known distances and from evolutionary tracks, and the results from 
both approaches are compared. 

Results. The major result of the present investigation is that the derived spectroscopic distances depend crucially on the applied reddening 
law. Assuming either standard reddening or values based on radio-H/3 extinctions, we find a mean distance of 9.0 ± 1.6 kpc and 12.2 ± 2.1 kpc, 
respectively. An "average extinction law" leads to a distance of 10.7 + 1.2 kpc, which is still considerably larger than the Galactic center distance 
of 8 kpc. In all cases, however, we find a remarkable internal agreement of the individual spectroscopic distances of our sample objects, within 
±10% to ±15% for the different reddening laws. 

Conclusions. Due to the uncertain reddening correction, the analysis presented here cannot yet be regarded as a consistency check for our 
method, and a rigorous test of the CSPN evolution theory becomes only possible if this problem has been solved. 

Key words. Stars: atmospheres - Stars: fundamental parameters - Stars: winds, outflows - Stars: distances 



1. Introduction 

Planetary nebulae (PN) in the direction of our Galactic bulge 
attract much interest, in particular because of their abundances. 
They are also important, however, because the Galactic bulge is 
one of the places in the universe with a reasonably well-known 
distance; for that reason it is a natural ground for testing meth- 
ods of PN distance determination. Our purpose in this paper 
is to test so-called spectroscopic distances, i.e., distances ob- 
tained from a determination of the basic atmospheric parame- 
ters of PN central stars (7^, log g) using model atmospheres to 
interpret the observed central star spectra. 

An alternative idea is to take good spectra of PN cen- 
ttal stars in the Magellanic Clouds which until now has only 
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been done in the far-UV faerald & Bianchil Eo04). To perform 
such a campaign in the optical has to await the advent of new 
technologies and instrumentation, e.g., adaptive optics at opti- 
cal wavelengths and/or space-born faint object high resolution 
spectroscopy. High specttal resolution is necessary, because 
unlike massive O stars, PN central star spectra are contami- 
nated by strong nebular emissions affecting the main H and He 
diagnostic absorption lines. In view of this, the Galactic bulge 
has the essential advantage over the Magellanic Clouds that it is 
much closer to us. It also has a disadvantage: the interstellar ex- 
tinction is much higher towards our bulge, and we do not fully 
understand the properties of the interstellar me dium in that di 



rection. We refer the reader to the study by IStasinska et al 



(1992), who compared Balmer decrement extinctions versus 
radio-HyS extinctions and concluded that for most distant PN 
the ratio R of total to selective extinction must be lower than 
the canonical value of 3. 1, and suggested a value of R = 2.7. 

However, for the moment the advantage of a smaller dis- 
tance more than compensates for the extinction disadvantage, 
and so we started a project made possible by the advent 
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of the high-resolution spectrograph HIRES at the Keck 10- 
m telescope (Mauna Kea, Hawaii). Spectra of several bulge 
central stars were taken from 1994 to 1996, and a prelimi- 
nary description of th e best of these spectra was presented by 



Kudrit zki et al.l (119971) . 



The 10-year delay between observations and the present pa- 
per can be attributed to the lack of confidence in the last genera- 
tion of model atmospheres. However, plane-parallel, metal-free 
non-LTE models are no longer the state of the art, and nowa- 
days we are able to produce much more sophisticated models 
that take into account blocking and blanketing by numerous 
metal lines, and feature hydrodynamically plausible winds in 
an extended, expanding atmosphere. The analysis of the objects 
performed in the present paper is based on this new generation 
of model atmospheres, and it turned out that the determination 
of stellar atmospheric parameters (e.g. T e g and logg) appears 
to be robust enough for spectroscopic distance determinations. 

Spectroscopic distances as determined here are not inde- 
pendent of stellar evolution theory; we can make the jump from 
T e a and log g to mass and distance because we assume that the 
relation between core mass and luminosity, and the post-AGB 
evolutionary tracks, are correct. In summary, what we are test- 
ing here is whether or not the combination of post-AGB evolu- 
tion theory and the presently best available model atmospheres 
can successfully predict the distance to the Galactic bulge from 
high-resolution optical spectra of a sample of five bulge PN 
central stars. We present a rigorous error analysis for this test. 

This paper is organized as follows: In Section|2]we comment on 
the selection criteria for our sample, the observations and data 
reduction. Section [3] comprises a brief summary of the NLTE 
atmosphere code used for the analysis, the method applied to 
derive stellar and wind parameters from the spectra and some 
comments on the profile fits and spectroscopic results for the 
individual objects of our sample. In Section H] we investigate 
those parameters which are not directly measurable from the 
spectra (particularly masses and radii), by introducing two dif- 
ferent analysis methods, based either on evolutionary tracks or 
adopted distances. We compare the results, with particular em- 
phasis on the involved error sources and the obtained degree of 
precision. We present the spectroscopic distances, and investi- 
gate the wind-momentum luminosity relation for our sample. 
In Section [5] we discuss the results and comment on several 
recent developmen ts like the new trigonometric parallaxes by 



Harris et al.l (12007). Section[6]summarizes our conclusions and 



provides future perspectives. 

2. Object selection and observations 

The PN central stars (CSPN) to be observed were selected 
among the brightest in the direction towards the Galactic cen- 
ter, according to the following criteria: 

(1) located within 10 degrees of the direction to the Galactic 
center; 

(2) angular size smaller than 10 arc sec, to reduce the probabil- 
ity of foreground PN; 

(3) only low-excitation PN, to ensure low T e R, because in that 
way we know the stars are evolving at high luminosities 



towards higher 7^; in addition, at low enough it is eas- 
ier to determine T e ff from the ionization equilibrium of He i 
and He n stellar features. 

For comparison with the previously studied "solar neighbor- 
hood" CSPN, we selected the central star of He 2-108. We 
have chosen this particular object because it has been ana- 
lyzed by different groups and, more importantly, using differ- 
ent methods. Previous analyses performed in a similar way as 
our approach result ed in a s pectroscopic distance of 6.0 kpc 
dMendez et all 19921 . see also lKudritzki et alJI 19971 and 2006), 
whereas alte rnative approache s were either consistent with 
these values ( Napi wotz kill2006l) or gave considerably different 
results (Pauldrach et al .1l2004 (see SectionE). 

The bulge CSPN spectrograms were obtained with the 
10 m Keck telescope, Mauna Kea, Ha waii, using the high- 



lUm JvecK telescope, Mauna Jvea, Hawaii, using 
resolution echelle spectrograph HIRES dVogt et al.l l 1 



1994). The 



spectra cover the range from 4250 A to 6750 A with a spec- 
tral resolution of 0. 1 A. Individual exposures were between 
100 s and 1800 s, and several exposures were made for each 
of the objects, cf. Table Q] In order to avoid excessive con- 
tamination of one order by neighboring ones, particularly at 
shorter wavelengths, we selected the rather short slit "C5" 
(7. Ox 1.148 arcsec) for our observations. This ensured a clearly 
defined inter-order minimum and reduced the effect of scattered 
light from strong nebular lines, which produced a halo affecting 
in some cases more than one echelle order. 

The CCD image contained 31 echelle orders, spread 
over 2048 x 1024 pixels. During each observing night dark 
frames, flat fields, and Th-Ar comparison spectra were ob- 
tained immediately before and /or after the object exposure. 
The spectrograms of the central star of He 2-108 were ob- 
tained with the 3.6 m telescope and Cassegrain echelle spec- 
trograph (CASPEC) of the European Southern Observatory, La 
Silla, Chile. For a de scription of the reductions we refer to 
Mend ez et al. (1988b). From these spectrograms we used HyS, 



He i Ai4471, 4387, 4922 and He n A4200. 

Additional spectra of the central star of He 2-108 were 
taken with the 3.5 m ESO New Technology Telescope (NTT) 
and EMMI Cassegrain spectrograph at the Cerro La Silla, 
Chile, in March 1994. From these spectrograms we used Hy, 
Ha, He i AA1 1 3 and He n and 454 1 . 

The usual calibration exposures (flat field, comparison 
spectra) were obtained. For brevity we will not explain the 
reduction of these simple spectra but will concentrate on the 
echelle spectra only. 

HIRES spectra reduction was performed with IRAFQ us- 
ing the following procedure: First we removed hot pixels and 
cosmic rays and subtracted the dark frames. Thereafter, the run 
of the order maxima was fitted, and a pixel shift between ob- 
ject and flat field frames was identified. Next, the object images 
as well as the comparison spectra were corrected for the nor- 
malized flat field images and the individual echelle orders of 



1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



Table 1. Summary of observations (Galactic bulge objects) 
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number of 


total exposure 


object 


observations 


time (seconds) 


He 2-260 


9 


6700 


M 1-37 


6 


12000 


M 1-38 


2 


3300 


M2-12 


5 


5500 


M2-33 


9 


8080 



the object spectra were extracted according to the fitting func- 
tions determined previously. Then the wavelength scale was 
calibrated using the comparison lines, which ensured an un- 
certainty of about 0.02 A. 

Finally, the extracted spectra were rectified. This process 
turned out to be quite delicate because of the high resolution of 
the spectra and the contamination by nebular light. Due to the 
high resolution, some of the broader stellar lines like Ha were 
spread out over more than one echelle order, which made the 
identification of the stellar continuum in this range rather dif- 
ficult. This challenge was met by interpolating the continuum 
run of neighboring echelle orders. With this method, we were 
able to rectify all echelle orders. 

The light of the nebula not only contaminates the stellar 
line profiles but also the stellar continuum. Therefore, an ad- 
equate correction is needed in order to obtain the real stellar 
line profiles, which (in case of absorption) become stronger in 
relation to the normalized continuum after accounting for the 
nebular continuum. When exposing the spectrograms, the neb- 
ulae were not fully covering the slit, and the insufficient spatial 
resolution of the nebulae prevented a simple subtraction. To 
correct for the nebular continuum passing through the slit, we 
ap plied the followin g procedure: 

Pottaschl d 19831 p. 84) showed that the ratio Ra of the inte- 



grated flux F'), eh in the H/3 line to the nebular continuum, F " eb , 



Rr = 



7 neb 



7neb 



const 



(1) 



only weakly varies as a function of the temperature and elec- 
tron density in the planetary nebula, as long as the consid- 
ered range is not too large. We adopt a representative value 
of Rp = 2000 A, as we expect the nebular electron tempera- 
tures to be relatively low, due to the low temperatures (28 000 K 
to 39 000 K) of our central stars. Since the spectrograms have 
already been normalized for stellar plus nebular continuum, 
the correction for nebular continuum must be performed us- 
ing the equivalent width Wp of the nebular H/3 line instead of 
the flux Fp eh . The corrected flux of the stellar continuum then 
follows as 



1 



Wp 
Rh 



(2) 



which, again, had to be normalized after the correction. 

A similar calculation can be made using the Ha and 
Hy emission lines, if we adopt the Balmer Decrement. For 
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Fig. 1. Famous emission lines at AA4485, 4504 in the spectrum 
of M 1-37 (see text). 



cool nebulae we h ave Hy/HyS = 0.47 and Ha/H/3 = 2.86 
(Ost erbrockl Il989l) . From this and the measured equivalent 
widths of the three emission lines, we found the average val- 
ues for the correction term, Wp/Rp (see Table©. 

2.1. Spectral morphology 

The central stars in our sample turned out to be all of the H- 
rich variety, i.e., with clearly visible H Balmer lines, implying 
a roughly normal He abundance, and describable in a first ap- 
proximation as late-O-type spectra, with the only exception of 
M 2-33, which has a somewhat hotter central st ar, Q5f(H) (for 
a general attempt at central star classification, see Mended 1991 
and references therein). 

In addition, it was immediately apparent that these stars 
show different degrees of wind contamination. He 2-260, 
M 2-12 and M 1-38 show O central star spectra dominated by 
stellar absorptions, while M 2-33 shows Hen /I4686 in emis- 
sion, traditionally classified as Of. 

Finally, we have the spectacular spectrum of M 1 -37, with 
rather weak He n /14541 in absorption and He n /14686 in emis- 
sion, but with impressive P Cygni profiles in Hei /L14471, 
4922, 4387, 4713, and especially in Sim AAA552, 4568, and 
4575 (see Figure|6j. The spectrum of this central star is unique 
and we are forced to classify it as peculiar. In fact, the Si m lines 
are so strong that one might be te mpted to classify this star as a 
low-gravity B0 or B2 (see, e.g., Walborn & Fitzpatrickl 199dh . 
which is of course impossible since He n features are present. 



Table 2. Average values of Wp/Rp, required to correct for neb- 
ular continuum contamination (see text). 



object 


He 2-260 


M 1-37 


M 1-38 


M2-12 


M2-33 




0.037 


0.020 


0.032 


0.030 


0.042 
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Strong Si m at higher temperatures is already a hint of a strong 
Si overabundance or at least of very specific, Si m enhancing 
atmospheric properties. 

Another interesting detail is the presence of very strong 
emission lines at /L14485, 4504 (FigureQ]). These famous lines, 
frequently seen in late-O supergiants and CSP Ns (see, e.g.. 
Swings & Strove! 1940l lWolff 1963, iHeadl 19771 and iThackeravl 



1977b . have been attributed to S iv bv lWerner & Rauchl d20oTb 



3. Determination of stellar and wind parameters 

3.1. The model atmosphere 

For the determination of the atmospheric and wind parame- 
ters of these stars we use the non-LTE, spherically symmet- 
ric model atmosphere code fastwind (Fast Analysis of STellar 
atmospheres with WINDs ), which was first introduced by 
Santolaya-Rey et alJ (I1997I) . Since then many improvements 
have been implemented, most notably the inclusion of line 
blocking/blanketing effects and the calculation of a consis- 
tent temperature structure by exploiting the condition of flux- 
conservation in the inner par t of the atmosphere and the elec- 
tron thermal balance (e.g., iKubat et al.l Q999) in the outer 
part of the atmosphere. An in-depth explanation of the meth- 
ods used in the code as well as a comparison of representa- 
tive results with those from alternative NLTE codes cmfgen 
dffillier & Millerll 19981) and WM -basic jPauldrach et alJbOOlh 



can be found in Puis et al 



(120051) . 

A number of spectroscopic investigations of CSPN have 
been performed u sing fastwind without blanketing (e.g., 
Kud ritzki et al.l 1997) and with blanketing and wind clumping 
( Ku dritzki et al.ll2006b . where the possibility of including wind 
clumping in the atmospheric treatment has been incorporated 
(and carefully tested) re cently. Our ap proach follows closely 
the original approach by ISchmutzl (1 1 995b . i.e., we assume that 
clumping is a matter of small-scale density inhomogeneities 
in the wincQ, which redistribute the matter into clumps of en- 
hanced density embedded in a rarefied, almost void medium. 
The amount of clumping is conveniently quantified by the so- 
called clumping factor, f c \ > 1, which is a measure of the over- 
density inside the clumps (compared to a smooth flow of iden- 
tical average mass-loss rate)H Remember that all processes that 
are controlled by opacities depending on the square of the den- 
sity (such as hydrogen recombination lines and bf/ff continua 
in hot stars) are preserved as long as the apparent mass-loss 
rate, M Jfa, remains constant, irrespective of the individual 
values of M and f c i themselves (at least as long as f c \ is spatially 
constant throughout the corresponding line/continuum forming 
flow). 

In the present analysis we considered the possibility of 
clumping for some of the objects (namely He 2-108, M 1-38, 
M 2-12 and M 2-33), in order to explain certain inconsisten- 



2 in contrast to large- scale inhomogeneiti e s suc h as co-rotating 
interaction zones (e.g., ICranmer & Owockil 1 19961 and references 
therein). 

3 An alternative description is based on the volume-filling factor, 
/v = /d 1 - 



cies between Ha and Hen 4686 arising in unclumped models 
(see Section l3~3T >. 

In a first approach, we used a rather simple stratification 
(but see IPuls et al.l i2006) for the clumping factor, allowing for 
an unclumped lower atmosphere (typically until r = 1 .01 7?*, 
where R* is located at Tr oss = 2/3), and a linear increase of 
the clumping factor until we reach its maximum (fit-) value 
(fd = 7. ..30) at r — 1.1 R„ (corresponding to typically 10% of 
Vco), assuming that the wind-line forming region is uniformly 
clumped. In this way, we account for the fact that any instabil- 
ity (in particular the line-driven one, which is generally con- 
sidered as the driver for wind clumping) needs a certain time 
to evo lve and to become non-linear (e.g 



Runacres & Owocki 



2002). Although there are observational indications (from UV- 



spectroscopy) that clumping might star t from the wind-base on 
(Hill ier et al.ll2003t iBouret et al.ll2.Q05b . so far we have found 
no problems with our simplified approach. 

3.2. Spectral analysis 

Projected equatorial rotational velocities v sin i (in the range 
from 50kms _I to 100 km s -1 ) were determined from the He- 
lines while fitting the theoretical profiles to the observed ones. 
Although rotation might not be the only effect for line broaden- 
ing3, the values derived in this way yield a fair estimate of the 
upper limits of v sin i. The ra/cro-turbulent velocity has been 
adopted as v m j c = lOkms -1 . Even though this parameter is not 
directly measurable (again, due to missing metallic lines), this 
value is consistent with evi dence from O/early B-stars in the 
considered parameter range dVillamariz & Herreroll2000b . 

The atmospheric parameters were determined from the hy- 
drogen Balmer lines. Hy, H/3, and Ha, as well as the Hei sin- 
glets AA4387, 4922, the He i triplets AA4471, 4713 and the He n 
,M4541 and 4686 lines (in case of He 2-108, Hen ,14200 was 
available as well). In a first step, we compared the observed 
profiles to synthe tic ones from our grid of OB star atmospheres 
( Puis et al.l2005 ), leading to rougrlj estimates for T^f, logg, th e 
wind strength parameter Q = M/(v«,S) 3/2 dPuls et al.lll996h . 
and the helium abundance Yu e = NHe/Nu- 

To improve the solution, all atmospheric parameters (in- 
cluding the clumping factor) have been fine-tuned to allow for 
the best possible fit, where the specific radius of the model 
was chosen from the evolutionary tracks (see Sec tion l4. 1.1 l and 
Figure [TTTl. updated in parallel with the progress of the line fits. 
The final model was re-run with a radius adjusted to the \ogg 
corrected for rotational velocity (see Section |4. 1 . ll for details). 

At first, the line profiles of Hy and H/3 were used to deter- 
mine log g for the possible range in T e ff as resulting from our 



4 E.g., macro-tu rbulence was fou nd to be significant at least in B- 
type supergiants dRvans et alj|2002h . but cannot be determined here 
due to missing weak metallic lines 

5 Note that the considerable differences in radii between normal 
OB-stars and CSPN (factors of 5 to 10) deteriorate any estimates re- 
lying on the usual scaling relations of wind diagnostics when using a 
grid of "normal" stars to investigate central stars: the wind conditions 
do not depend on the wind-strength parameter alone, but also on quan- 
tities with a dif ferent scaling behavior, such as the density (see, e.g., 
IPuls et alj|2005t) . 
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Fig.2. M 2-33: Hy at different gravities. r eff = 39 000, T He = 
0.1, v sin/ = 80 and wind strength as derived from our com- 
plete analysis (see Table [3j. Note the difficulty in deciding on 
the best fit. 



crude estimates. The effect of different gravities on Hy is dis- 
played in Figure|2] Then, the effective temperature was derived, 
from the ratio of He 1 to He 11, whereas the absolute strengths 
of the helium lines define the helium abundance, Fne- Finally, 
the wind strength parameter, log Q, and the clumping factor f c \ 
were determined, using all profiles, in particular Ha and He 11 
4686. 

In case of strong mass-loss (well-developed Ha), v M could 
be determined from the width of Ha (we used v^/v^ ~ 2 to 3 
as a starting point), whereas in cases of weaker M the synthetic 
profile reacts too weakly on variations of v M to allow for such 
a procedure. For these stars Vco was interpolated from the Table 
provided bv lKudritzki et alj (11997b , who used terminal veloci- 
ties for CSPN of roughly similar parameters derived from UV 
resonance line profiles. 

For all objects, we also checked the available silicon lines 
for consistency. The c orresponding model atom is the same as 
used and described bv lTrundle et alj (|2004) in their analysis of 
SMC B supergiants. 



3.3. Wind clumping 

During our first analysis of the spectra, for two of our objects, 
M 1-38 and M2-12, it was impossible to fit Ha and Hen 
/14686 at the same mass-loss rate simultaneously (the predicted 
He n emission was by far too large, compared to the observed 
absorption profiles, see Figure[8] dotted profiles). 

Since at that time evidence had a ccumulated that WR and 
OB-star winds might be clumped (cf. IPuls et al.ll2006l and ref- 
erences therein), and since it was to be expected that clumping 
should affect Ha and He 11 in a different way for cooler objects, 
we investigated to what extent this process is able to solve the 



discussed problem. The assumption of a clumped wind leads to 
a satisfying solution. 

Independent of this work, the same problem and solution 
was found in the re-analysis of a different sample of CSPN, and 
a detailed discussion of the involved physics can b e found in the 
corresponding publication (IKudritzki et al . 2006). In brief, and 
as outlined above, small-scale clumping preserves the profiles 
of p 2 -dependent lines, as long as the apparent mass-loss rate, 
M Jf^, remains constant. Though for almost the complete OB 
star range Ha follows this scaling (since H 1 is a trace ion), He n 
becomes the dominant ion (or close to this) for late O-type ob- 
jects, and the corresponding opacities scale (almost) linearly 
with density. In such a case then, the corresponding line pro- 
files (in particular, He 11 /14686) do not react on wind-clumping, 
but depend on the actual mass-loss rate, M, alone. Thus, a 
clumped wind with a mass-loss rate lower than derived from 
unclumped models can explain that He 11 /I4686 is still in ab- 
sorption (because the actual mass-loss rate is low), whereas the 
strong emission in Ha is a consequence of clumping. Figure[8] 
illustrates this behavior, where the solid grey profiles have been 
calculated with a mass-loss rate reduced by a factor of y^i 
(compared to the dotted solution), and a clumped wind with 
/ c i = 30 has been assumed. Note that all other lines analyzed 
here react only weakly to the presence/absence of clumping, 
either because they are formed very close to the photosphere, 
or because they are also dependent on p 2 (such as He 1 at cooler 
temperatures). 

Two additional comments: (i) For hotter objects 
(r e ff 35...38kK, depending on gravity), Hem becomes 
the dominant ion, and both Ha and He 11 /14686 depend on 
the apparent mass-loss rate. In such a case, the comparison 
of both lines does not allow us to conclude on the degree 
of clumping, (ii) If the observed Ha and Hen /14686 lines 
show a similar degree of emission at cooler temperatures, this 
strongly suggests that the corresponding wind is unclumped, 
i.e., /ci = 1. 

3.4. Error estimates for quantities derived from the 
spectra alone 

Thanks to the high quality of our spectra, fitting errors due to 
resolution limitations or instrumental noise do not play a role in 
our analysis. The major problem encountered here is the strong 
contamination of the hydrogen and He 1 lines by nebular emis- 
sion in their cores and beyond. In certain cases (particularly in 
Ha, e.g. M 1-37 and M 1-38) the contamination is so severe 
that the uncertainties in the derived parameters depend mainly 
on our inability to identify the stellar profile unambiguously. 

Apart from this principal problem, the second major source 
of errors is due to our "eye-fit" procedure (contrasted to au- 



tomated methods, e.g., Mo kiem et al.ll2005l) . which in certain 



cases might not be able to find the global optimum in parameter 
space. Because of the severe nebular contamination, however, 
an automated method is very difficult to apply, due to the prob- 
lems of defining the "true" stellar spectrum required for such a 
procedure. 



6 



P. Hultzsch et al.: CSPN in the Galactic bulge 



In the following, we will summarize the typical errors 
resulting from our spectral analysis. Errors of quantities re- 
lated to other diagnostics (masses, radii, distances and absolute 
mass-loss rates) will be discussed in Section[4] 

3.4.1 . Effective temperatures 

Depending on the quality of the helium line fits and the neb- 
ula contamination, the typical formal errors in T e ff are about 
+ 1 500 K. This uncertainty becomes larger in the region around 
r e ff = 32 000. . .36 000 K, if we account for the (known) prob- 
lem of a possible inconsistency between the predicted line- 
strengths of He i triplets and singlets (due to subtle line over- 
lap effects regardin g the Hei resonance line in the FUV, see 
Najarro et al.l l2006). In this tempera ture range, the triplet lines 
are the more reliable ones (see also iMokiem et al.ll20Q5h . and 
are consequently used as a main indicator for the determination 
of T eB . 



3.4.2. Gravities 

Because of the intrinsic coupling of log g and T e ff when deriv- 
ing the gravity from the Balmer line wings (typically, a change 
by 1000K corresponds to change in \ogg by 0.1, in the same 
direction), the precision of log g (which is of major importance 
in the present investigation) depends strongly on the precision 
of Teg (from the He ionization balance) which in turn requires 
a reasonable determination of the helium abundance, Yue- With 
respect to the available spectra and their nebular contamina- 
tion, this may not be warranted in all cases. Particularly for 
He 2-108, the uncertainties were large enough to persuade us 
to provide two different solutions for this object, at lower and 
higher T e s (for details, see Section l3.5.6l >. 

In most other cases, we adopt, in accordance with sim- 
ilar investigations, a typical error of Alogg = +0.1, since 
for larger differences the profile shape changes significantly. 
If the star is close to the Eddington limit, the effective grav- 
ity becomes smaller and smaller because of the increased ra- 
diation pressure (not only from Thomson scattering, but also 
from photospheric line pressure). The "true" gravity is given 
by gg ra v = geff + grad (excluding centrifugal forces for the mo- 
ment, but see Section !?. 1. Il l, where the Balmer line wings actu- 
ally "measure" g e ff because of their dependence on the electron 
density. As soon as g ral j ^ g e g, any uncertainty in the calcu- 
lated radiation pressure (intrinsic to the atmospheric models 
and strongly dependent on the completeness of the line list) 
results in a larger error for g erav , which we adopt in these cases 
byAlogg = 0.15. 

3.4.3. Wind parameters 

As outlined above, terminal velocities were derived, where 
possible, from the width of Ha. Particularly, we found Vco = 
500 ± lOOkm/s for M 1-37, and v m = 800 + 160km/s for 
M 2-33. Strictly speaking, both values can be considered as 
lower limits only. For the other objects we adopted Voo from in- 



terpolating between " known" terminal veloc ities measured for 
different samples (cf. iKudritzki et al.l Il997l) . and allowed for 
typical errors of 20%. 

The precision of the wind-strength parameter, log Q, on the 
other hand, depends on the validity of the derived velocity ex- 
ponents, p ( Puis et al. 19961) . To account for this fact, we have 
fitted log Q for different values of f3 + A/3, with A/3 = 0. 1 . In this 
case then, A log Q i s ty pically not larger tha n +0.1 (see also 
Markova et a"i]|2004l and lRepolust et alj|2004 . 



3.5. Comments on the individual objects 

In the following, we compare our best fitting synthetic spectra 
with the observed ones, for all H and He lines which have been 
analyzed. For two of our objects, He 2-260 and M 1-37, the 
Si line "fits" are displayed as well. Whereas for all but one 
of our objects the Si lines matched nicely the observations at 
solar abundance (with He 2-260 being a prototypical object in 
this respect), we found a very large Si-abundance for M 1-37, 
which is still lacking any theoretical explanation. 

Note that all of our objects display a rather strong Ha 
emission, immediately pointing to large wind-strengths (cf. 
Table O, which indeed are a factor of about 5 to 10 larger than 
for typical O and B supergiants. Further comments on this be- 
havior are given in Section l4~3l 

All results from our analysis including corresponding er- 
rors have been summarized in Table [3] 

3.5.1. He 2-260 

For this star, we obtained a convincing fit, displayed in 
Figures [3] and |U The observed and synthetic Hy lines agree 
well, i.e., logg is well determined within the nominal error. 
The second diagnostics for log g, H/3, would require a slightly 
larger value, but this would deteriorate not only Hy but also 
He i /U4387 and 4922 as well as Si m. 

Compared to the other objects, Ha shows the lowest emis- 
sion, which immediately points to a "low" mass-loss rate. Since 
He ii /i4686 is closely reproduced in parallel (no clumping re- 
quired to keep this line in absorption at the derived effective 
temperature of 28 000 K), and all other lines fit almost perfectly 
(including silicon, where Si in is matched at solar abundance, 
and Si n and Si iv are absent as predicted) we feel rather confi- 
dent about the stellar and wind parameters of our model. 

3.5.2. M 1-37 

Although not classified as a Wolf-Rayet star, this star shows 
well developed P-Cygni profiles in He i (Figure |5]l and Si m 
(Figure [6]), immediately pointing to an extremely dense wind 
(the densest in our sample). 

To fit the He i lines, an enhanced helium abundance of at 
least Yu e = 1 .0 is required, where even larger values cannot be 
excluded due to saturation effects. For lower helium content, 
on the other hand, a fit to the He i lines would require a larger 
mass loss rate, leading to too strong emission in Ha, H/3 and 



i.e., not accounting for any errors intrinsic to the synthetic profiles. Hy. 
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He 2-260: H and He i/He n lines analyzed in this study, together with the predictions from our best fitting model (grey). 



Although there is some (stellar?) absorption visible in the 
blue wings of Hy and H/3, we were not able to match these 
features, thus leaving the gravity rather unconstrained: Models 
calculated with different log g (and all other parameters kept 
constant) display only marginal differences in the spectra, 
which is not surprising since for such large wind-densities all 
lines are formed in the wind. However, if we changed not only 
log g but also R, (and M to conserve Q) according to the core 
mass-luminosity relation, differences becam e visible, becaus e 
not all quantities scale exactly with Q (cf. IPuls et al.1 120051) . 
The differences for A log g = +0.3 can be seen in Figs.|5]and[6] 
The largest reaction is found for He n /14686, where the emis- 
sion becomes stronger with lower \ogg and larger R„. Note, 
however, that this effect can be compensated by introducing a 
certain degree of clumping (see below). 

As we were not able to derive \ogg spectroscopically, we 
adopted \ogg ^3.4 from distance considerations (assuming 
that the star is located at w 8 kpc). Lower values for log g would 
place the star further away, higher values to closer distances. 



Moreover, we adopt quite a large error of Alogg = +0.3, 
thereby allowing for distances in between 5.3 kpc and 1 1.5 kpc 
(see Table This seems to be a reasonable range, compar- 
ing with our most distant object, M 2-12, located at 11. 4 kpc 
according to our analysis. 

The determination of T e s from the ratio of He i to He n pro- 
vides similar problems as the definition of the other parame- 
ters. The adopted value of 31 000 K is an upper limit, depend- 
ing alone on our fit to Hen /14541. A lower limit on T e ff can 
be derived from the absence of stellar Si n lines, which have 
been found to appear below T e ff = 29 000 K. To fit the Si in 
emission lines, an extremely large Si abundance is required (28 
times the solar abundance). According to stellar evolution mod- 
els th e Si abundance should n ot change, even in H-deficient 
stars (IWerner & Herwig 2006), and we have no explanation 
for this finding. Leuenhagen & Hamann ( 1998) claim solar Si 
abundance for similar objects but mention overabundances (up 
to 40 times solar) for some [WCL] stars. 
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Fig. 4. He 2-260: important Si lines present in our data set. Overplotted are the theoretical predictions for identical parameters 
as used in Figure[3](solar Si abundance). 



Both Ha and He n /14686 are consistent with an unclumped He n emission and require a higher T e ff to recover the observa- 

wind of the quoted strength, although the complete Ha line is tions, which is inconsistent with the other lines. Nevertheless, 

hidden within the strong nebular emission. Note that there is no (moderate) clumping cannot be excluded, since a lower grav- 

indication of clumping at the derived/adopted parameters of our ity (compensating for the decreased emission) is still possible, 
final model, since the presence of clumping would decrease the 
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Fig. 5. As Figure[3] but for M 1-37. The dotted profiles refer to models with different gravities, A log g = ±0.3 (stronger emission 
in He n/14686 for decreased gravity, see text). 



Note also that a clumped wind might be able to decrease the 
derived Si abundance. 



3.5.3. M 1-38 

All parameters can be derived with a reliable precision within 
the nominal errors, log g is well defined from the wings of Hy 
and H/3, and T e s and (solar) Yn e are equally well constrained 
from Hen AA454 1/4686 and He i ,14471. In contrast, however, 
we could not obtain a satisfying fit for the He i /1/14387/4922 
singlets (overestimated strengths, presumably due to the "sin- 
glet problem" outlined above), and we have discarded them 
from our analysis. 

This object is the first where we needed a clumped wind, in 
order to fit Ha and He n /14686 simultaneously at the derived 
effective temperature. Note that all other lines do not react on 
clumping, but are consistent with the derived "apparent" wind- 



strength, M ^Jfc\■ This value and the individual quantities for M 
and / c i « 20 were determined by a compromise to fit the wind- 
induced asymmetries seen in most lines, the emission profile in 
Ha and the absorption line profile of He n /14686. 

3.5.4. M 2-12 

Both Hy and H/3 are dominated by the strong nebular lines, and 
only the extreme wings become visible. This leads to a higher 
uncertainty in logg, with an error of Alogg = +0.2/ - 0.1. 

As above, the synthetic He i singlets might appear as too 
strong, but in this case the nebular contamination forbids a reli- 
able estimate on whether this is really the case. From the triplet 
lines and He n /1/1454 1/4686 then, we derive a well defined T e ff 
and helium abundance Yu e , within the nominal errors. We pre- 
fer a somewhat reduced helium abundance, Yu s = 0.08, in or- 
der to fit the globally weak He-lines. 
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Fig. 6. As FigureHl but for M 1 -37. A silicon abundance of roughly 30 times the solar one has to be adopted in order to reproduce 
the observations. Again, the dotted profiles refer to models with different gravities, A log g = ±0.3. 



Again, a clumping factor of f c \ m 30 is required to keep 3.5.5. M 2—33 
He ii /14686 in absorption. Further details on this problem have 

been outlined in Section [3~3l ^ or tms °bj ect > tne quality of both the spectra and the fits is 

satisfactory, and nominal errors apply for all parameters except 
for T e ff and thus logg (see Figure [9]). Since most of the Hei 
absorption lines are very weak and thus dominated by nebular 
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emission, the ratio of He i to He n cannot be determined well 
enough to ascertain the effective temperature precisely. We es- 
timate the corresponding error as Ar e ff = +2 000 K. 

Anyway, in this temperature range (roughly 39 000 K) He n 
is no longer a dominant ion, and the corresponding opacities 
scale with p 2 . According to our arguments given in Section [3~3l 
a unique determination of the clumping factor is no longer pos- 
sible, and we adopt f c \ = 1 . A comparison between such an 
unclumped and a clumped model is given in Figure [9] the dot- 
ted profiles refer to a clumped model with f c \ = 7, and the 
small difference in Heii/14686 result from the fact that Hem 
recombines in the outmost wind (v/vk, > 0.8). 



3.5.6. He 2-108 

This object has been previously analyzed in the literature, with 
two differing values for (and related quantities). A low tem- 
perature solution with T e ff = 35 000 K/ 34 000 K was deduced 



by Kudritzki et al. (1997/2006), from the optical, whereas a 
high temperature solut ion with T e ff = 39 000 K was found by 



Pauldrach et al 



d2004 . on the basis of the Fem/iv ionization 
equilibrium derived from the UV. 

Given the low quality of the available spectra, none of these 
solutions can be excluded. Three of the four available He i lines 
appear to be weak, thus supporting the high temperature solu- 
tion, whereas He i /14471 seems to be much stronger, thus favor- 
ing the lower T e g. Rather than trying to find a compromise, we 
decided to show both solutions, and to follow the consequences 
throughout the remainder of this paper. 

Our low temperature solution (Figure \TU[ solid) has been 
optimized for He i AAA1 1 , and results in T e s = 34 000 K, identi- 



cal to the value derived by Kudritz ki et al. (2006). In this case, 
both He i singlets ,L143 87/4922 and the Hei ,14713 triplet are 
overestimated, whereas for the high temperature solution (dot- 
ted) it is basically the other way round. Note that the singlets 
at 7rff = 34 000 K might be overestimated, due to the singlet 
problem. At least at T e # = 39 000 K, however, this problem is 
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Fig. 8. As Figure [3] but for M 2-12. The dotted profiles display the predictions from an unclumped model at lar; 
rate. Note that the only difference regards He n ,14686 (see Section l3~3l . 



4927. 4932. 

rer mass-loss 



no longer relevant, and the agreement between observed and 
predicted singlets cannot be disputed. 

The errors for our low temperature solution (from Hei 
,14471) have been estimated as +1 000/-2000K for T eS , and 
a standard error of Alogg = 0.1. In this case then, the helium 
abundance must be slightly enhanced, Yn e = 0.15. 

From the simultaneous analysis of the Balmer lines and 
Hen /i4686, we found an interesting result on the clump- 
ing properties of the wind. In contrast to the other objects, 
He 2-108 seems to require larger clumping factors in the in- 
ner wind, compared to the outer one. First, note that, at T e ff <: 
34 000 K, there is no discrepancy between He and He n /14686 
(in agreement with the discussion given in Section [3~3l ). In con- 
trast, however, there is now an inconsistency between Ha and 
Hy/H/3. 

For the low temperature solution, we have "normalized" the 
(unclumped) mass-loss rate to fit Hor/Hen ,14686, and obvi- 
ously the cores of Hy and H/3 are predicted as too deep. Thus, 



clumping must be larger in the inner wind, to allow for a refill- 
ing of both lines due to increased emission. 



For the high temperature solution, an unclumped model de- 
signed to fit Ha predicts too much He n ,14686 emission (this 
inconsistency cannot be cured by clumped models, because of 
the high temperature), and, again, too little emission in Hy and 



Note that He 2-108 is the only object in our sample which 
indicates such a stratified clumping. Because in the present 
analysis we are considering only depth independent clumping 
(see Section [XTl i, we provide (both for the high and the low 
temperature solution) a clumping factor of unity in Table [3] 
consistent with the wind-lines, but are aware of the fact that the 
real mass-loss rate(s) might be smaller than indicated. 
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Fig. 9. As Figure[3] but for M 2-33. The dotted profiles refer to a clumped model with f c \ = 7 (see text). 




4. Deduced Parameters: Masses, radii, distances 
and wind-momentum rates 

4.1. Two different methods 

Having presented those results that can be derived from spec- 
troscopy alone, we will concentrate now on the "remaining" 
parameters that are not directly deducible from our analysis, 
since they depend on the stellar radius, R*. Two different, inde- 
pendent approaches will be applied to obtain these parameters: 
"Method 1" uses evolutionary tracks for central stars to provide 
masses, which, together with log g, result in stellar radii and fi- 
nally spectroscopic distances. "Method 2", on the other hand, 
uses adopted (bulge-)distances to calculate radii and depending 
quantities. We will give a detailed error analysis and compare 
the results and inherent uncertainties of both methods. This will 
in particular clarify (at least for our sample of objects) the ques- 
tion raised in Section [1] namely in how far the combination 
of post-AGB evolution theory and state-of-the-art model atmo- 



spheres can successfully predict distances from high-resolution 
spectroscopy. 

So far, the spectroscopically determined gravities (from the 
Balmer-line wings, denoted by logg spec in TableO are effective 
values, since centrifugal forces lead to an outward accelera- 
tion of the atmospheric layers, thus reducing the "true" gravity, 
log g true . In order to obtain this quantity (which has to be used in 
our subsequent analysis), we have to correct log g spec by adding 
the centrifugal acceleration. This has been accomplished in the 
"usual" way by means of the derived projected rotational ve- 
locities, v sin L and the average v alue of sin i (for details, see 
Appendix in iRepolust et al.ll2004l and references therein). As 
already outlined in Section 13.21 due to missing metallic lines 
the macro-turbulent velocity, v maC ro, was not measurable for our 
sample, which means that the quoted values of v sin i and thus 
the true gravities have to be considered as upper limits. Note, 
however, that despite the rather low rotational speeds the dif- 
ferences between log g spec and log g tme are quite significant (of 
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Fig. 10. He 2-108, H and He line fits. Note the lower quality of the observed spectra, compared to the other objects. Solid: low 
r e ff solution, for an unclumped model and a mass-loss rate required to fit Ha and Hen /14686. Dotted: high T e ff solution, again 
for an unclumped model designed to match Ha. (see text). 

the order of 0.05 dex), which is due to the small radii of our 4.1 .1 . Method 1 
objects (centrifugal acceleration oc v^ ot /R*). 



Here we appl y the method to derive spectroscopic distances as 
described by iMendez et al. I dl988bh . We use the evolutionary 
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Table 3. Results of our analysis for the complete sample, using either "Method 1" or "Method 2" (see text). logg spec refers to 
the gravity as derived from spectroscopy, and log g tme is an estimate (upper limit) for the true gravity, corrected for centrifugal 
acceleration and used to calculate related physical properties. Wind strength parameter log Q calculated with M in units M /yr, v«, 
in km/s and R t in R . Modified wind-momentum rate, D mom , calculated in cgs. The different distances resulting from "Method 1" 
refer to the different extinctions as given at the top of the table. All other parameters which are dependent on extinction refer to 
the "mean" value, Ay me an (see text). 
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tracks in the log r e ff-log g diagra for CSPN from Bloecker 



7 Actually, those tracks are tabulated as a function of luminosity and 
effective temperature, which can be easily converted into the required 
log g vs. 7J. ft diagram due to given masses. 



(1 1995b and ISchonbernerj d 1 989l) (Figure HB, which defines the 
stellar mass, M», by interpolating the spectroscopically deter- 
mined r e ff and log g tme values between the evolutionary tracks 
of different masses. This method is feasible for CSPN that are 
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of 25% in the distance, which accounts for th e depth of the 



4.60 4.55 4.50 

og Teff 



A.A5 



Fig. 11. The positions of the CSPN in the log T e a - \ogg dia- 
gram, with error bars according to the derived uncertainties, to- 
gether with theoretical post- AGB evolutionar y tr acks for stars 
of five different masses (from [Bloeckerlll995l and Schonberner 
19891 masses in units of M ). L, H, K and P denote vari- 
ous results for He 2-108: our tow and high temperature so- 
lution, the analysis by /Tudritzki et al. (2006) and by Pauldrach 
et al. (2004). Remember that the gravity of M 1-37 has been 
adopted. 



in the phase of almost constant luminosity, because the evolu- 
tionary tracks for different masses are unique in this regime. 

Since the centrifugal correction term involves the stellar ra- 
dius, an iterative scheme is required: starting with some guess 
for R*, we obtain M», a new radius via logg tvue , which is up- 
dated subsequently and a new mass derived. The iteration pro- 
cess is continued until mass and radius have converged. In 
Table [3] we provide these converged quantities and the corre- 
sponding true gravities. Note that the masses and radii derived 
in this way depend on the validity of the tracks (and on the ac- 
curacy of the spectroscopic analysis), but not on the extinction. 

4.1.2. Method 2 

Here we exploit the fact that the observed objects are 
supposed to be located in the Galactic bulge, follow- 
ing our selection criteria (Section [2]). We adopt the well- 
defined distance to the Galactic center , Rq, of roughly 8 kpc 
(7.94 + 0.42 kpc. Eisenhauer et al.l2003l) for all our objects, ex- 
cept for He 2-108. 

For the latter one, we adopt a distanc e of 4.6 kpc which 



is consis tent with the value p rovided by Men dez et al.l (|1992, 
and also iKudritzki et al. 2006h . 5.8 ... 6 kpc, if we account for 
"our" extinction (Ay mean , see below) and visual magnitude V, 
together with their stellar parameters. Admittedly, this distance 
is far from being as well determined as Rq, and moreover it 
has been determined by using "Method 1". In so far, the corre- 
sponding results should be considered as a comparison rather 
than a test. 

To check the error propagation of "Method 2" and to com- 
pare it with "Method 1", in all cases we adopt a typical error 



Galactic bulge and follows the considerations by Mendez et al 
dl988bh . 

From the distances, the visual magnitude, V (cf. Table [3j 
all from Tvlen da et al. Il99ll) and the extinction, Ay (see be- 
low), we find the absolute magnitude, My, which together with 
the theoretical V-band fluxes (convolved with the correspond- 
ing filter transmission function) allows us to derive the stellar 



radius, R„, fol lowing Kudritzki (1980, see also Eqs. (1,2) in 
Repolust et al.ll2.004b . The errors inherent to all related quan- 



tities (mass, luminosity etc.) are thus dominated by the uncer- 
tainty in distance and extinction, except for M 2-33, where the 
large uncertainty in V becomes decisive. 

4.1 .3. Extinction 



Stasinska etail dl992h 



con- 



As outlined in the introduction, 
eluded, from a comparison of Balmer decrement extinctions 
versus radio-H/3 extinctions, that for most distant PN the cor- 
responding line of sight extinction law should be significantly 
different from the standard law (for the diffuse ISM). They ar- 
gued that the ratio Ry of total to selective extinction (in the 
V-band) must be lower than the canonical value of 3.1, of the 
order of Ry = 2.7. Similar conclusions were reached more re- 



cently bv lUdalskil d2003l) and lRuffle etai] d2004 . 

Using previous results by ISeatonl dl979l) . the standard vi- 
sual extinction can be related to the logarithmic extinction at 
Hy6 determined from the Balmer decrement, C opt , via 



Ay optical = 2.16 C opt ("standard" reddening). 



(3) 



The measurements of C ra d (the logarithmic extinction at H/3 de- 
termined from the radio fluxes) for a large sample of PN by 
Stasinska et al., on th e other hand, imply Ry = 2.63 (using the 
relations provided bv lNandv et alJll975l) and finally 



A v radio = 2.10 C la d (lower reddening for distant PN). (4) 
For all our objects, C op t > C rar i (values taken from Stasinska et 



al. and Tvlenda et al. 1992), and we have used both the stan- 



dard reddening to obtain A^opticai and the lower value Ay radio 
in our further investigations to estimate a maximum and min- 
imum effect. Moreover, we have used also a linear mean, 
Aymean = 0.5 (Ay opt i ca i + Ay ra dioX with larger errors to account 
for the present uncertainties. In Table [3] we display the differ- 
ent distances determined by "Method 1" resulting from these 
three values of Ay, whereas, for brevity, we display only the 
solutions relying on Ay mean for "Method 2". 



4.2. Comparison of Method 1 and 2: stellar 
parameters and distances 

Regarding the error propagation of the quantities observed or 
determined directly from the spectra (see first part of Table [3} 
into deduced physical quantities such as radius, mass, luminos- 
ity and mass-loss dependent quantitie s we refer the read er to 
the detailed di scussion provided bv lMarkova et alj ([2004) and 
Repo lust et ail d2004 . 



P. Hultzsch et al.: CSPN in the Galactic bulge 



17 




0.6 0.8 1. 
mass (Method 1 ) 




26 27 
og D mom (Method 1 ) 



•ft i /M 1 -37 



— < 


> 










-< 


> 




< 


>— 





2 3 4 5 

radius (Method 1 ) 




5 10 
distance (kpc, Method 1) 



' 5 



Fig. 12. Comparison of various parameters of our CSPN sample, derived by "Method 1" (based on evolutionary tracks) and 
"Method 2" (based on adopted distances), using the "mean" extinction, Ay mean (see text): masses in units of M , radii in R Q , 
modified wind-momenta in cgs and distances in kpc. L and H refer to the low and high temperature solution for He 2-108, 
respectively. The solid line displays the one-to-one relation. Note that most of the results agree within the error bars, but that the 
errors inherent to "Method 1" are considerably lower than those relying on "Method 2". Since the gravity of M 1-37 has been 
adopted from distance arguments, all parameters agree by definition. 



4.2.1 . Mass and radius related quantities 

Figure [12] (upper left panel) shows a comparison of the de- 
rived masses from the two methods. Obviously, the error 
bars from "Method 1" are much smaller than those resulting 
from "Method 2", which of course is related to the fact that 
"Method 2" involves the rather large uncertainties in /?» (due 
to uncertainties in distance and extinction), but, even more, 
because M, depends quadratically on radius. In contrast, for 
masses based on evolutionary tracks the dependence on errors 
in log g and T e s should be minor since the slope of the tracks is 
fairly parallel to the r e ff-log g iso-contours of equivalent width 
for the Balmer lines formed in the photosphere. 

Relying on this notion, we defined the mass uncertainty re- 
sulting from the adopted errors in log g and TJg in the following 
way: We determine M» for the pair (log g + A log g, T e ff ) as well 
as for the pair (log g, T e s + ATkf), and use the higher of both 
masses as our upper limit for M m . The lower limit is defined 
analogously, with opposite signs. We do not use combinations 



of (+Alogg, -Ar e ff) and vice versa, since such combinations 
can be excluded from the results of our analysis (remember that 
the errors in T e s and logg are not independent). 

The resulting uncertainties for M, are amazingly small; 
even the object with very large uncertainties in logg, M 1-37, 
has a fairly well determined mass, which is also related to the 
fact that for lower masses the dependence of M» on gravity 
is rather low. The mass range within our sample extends from 
0.58 to 0.78 M Q , with errors of about 15%. 

From these results, however, one might also argue that our 
sample is not completely representative. This is best seen from 
the fact that the mean stellar mass of our sample is about 
0.68M Q , which is considerab ly larger than the expected O.6M 
(see, e.g jLiebert et al ]l2005h . Such a discrepancy might be due 
to selection effects, because for our analysis we choose only ob- 
jects with a high enough signal to noise ratio. With such a bias 
towards more luminous objects, we might have selected also 
the more massive ones. 
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In contrast, using "Method 2", the masses are determined 
from the radius which enters quadratically and is dominated 
by the large uncertainties in distance and extinction or visual 
magnitude. If we consider the corresponding errors as already 
discussed and summarized in Table[3] we find errors in between 
30% to 45% for R, and uncertainties of the order of a factor of 
two (larger or smaller) with respect to M, . 

Comparing the luminosities shows a similar result. Except 
for M 1 -37, the errors in log L„ derived by "Method 2" are 
twice as large as those from "Method 1", whereas in the case 
of M 1-37 the errors are similar, due to the extremely large 
uncertainties in log g affecting the corresponding stellar radius 
derived by "Method 1". 

Concentrating on the parameters themselves (mass and ra- 
dius), Figure [12] shows that for most of the objects the derived 
values agree within the error bars, but that the specific val- 
ues can be quite different. Notably, the masses of three bulge- 
objects as derived from "Method 2" would be extremely low, 
whereas both solutions for He 2-108 result in implausibly high 
masses. Again, this is the outcome of the quadratic influence of 
a highly insecure stellar radius. 

4.2.2. Distances 

Provided that our atmospheric models do not contain intrinsic 
errors and that the evolutionary tracks of the central stars are 
generally applicable to all objects of this group, the uncertain- 
ties in the spectroscopic distances are dominated by the errors 
in the extinction towards the central stars (plus the error in V 
for M 2-33). Remember that the derived masses are very accu- 
rate in this case, and that the errors in stellar radius result from 
a fairly small error propagation since R, oc yjM t /g. In order to 
account for the unclear situation regarding reddening, we as- 
sumed rather large errors for Ay optical (standard reddening) and 
^Viadio (Rv ~ 2.6), namely the mean difference between both 
quantities for our sample (AAy = +0.6 mag), whereas for the 
"mean" extinction, Ay mean , we assumed, because of the applied 
compromise, individual errors of AAy mean = +0.5 (Ay rat ji - 
Ayopticai) which are roughly a factor of two smaller. 

The resulting distances for Ay mean are displayed in 
Figure [T2] (lower right panel), and those using either Ay opt i ca i 
or Ay ra dio are tabulated in Table [3] At first note that for the 
supposed bulge objects we derive distances in the range of 
7.8 ...11.5kpc when using the "mean" extinction, whereas 
for the alternative solutions Ay opt i ca i and Ay ra di we derive dis- 
tances that are smaller and larger, respectively. 

Regarding the specific values, the distances derived by us- 
ing Ayopticai are best compatible with a bulge population. This 
becomes immediately clear if we assume that our sample is 
statistically distributed, and we calculate an (error-weighted) 
mean distance (leaving out M 1-37 since it has been placed 
inside the bulge, see Section l3.5.2l i. In this case, we obtain the 
values found in Table |4] Further consequences of these results 
are discussed in Section[5] 

Both distances that we obtain for He 2-108 (low and high 
temperature solution) are fairly similar, a nd the differences with 
respect to the re-scaled distances from iMendez et al. (1992) 



Table 4. Derived mean distances using different Ay disregard- 
ing object M 1-37 (see text). 



extinction 


^Vmean 


A V optical 


A y radio 


mean distance (kpc) 


10.7 + 1.2 


9.0 ± 1.6 


12.2 + 2.1 



and Kudritzki et al. ( 20061) (compare with the entry according 
to Ay mean ) are due to the different gravity resulting from our 
analysis and the applied centrifugal correction, which has been 
neglected previously. For this object, the influence of different 
extinction/reddening-laws is not as large as for the bulge ob- 
jects, because the extinction itself is low. 

4.3. Mass loss rates and wind momenta 

Since the diagnostics used (Ha and He n /I4686) constrain the 
wind-strength parameter, Q (or, even worse, only the product 
Q -\ffci), aiso the mass-loss rates and the (modified) wind mo- 
mentum rates depend on radius (and via M — Q (R„ Voo) 3 ^ 2 
and D mom = Q R 2 t Vol 2 . Figure [12] (lower left panel) compares 
the latter quantities as resulting from "Method 1" and "2". At 
least on the logarithmic scale, they are quite similar, though the 
errors from "Method 2" are larger, of course. 

Figure Q~3] displays our results together with those de- 
rived from other investigations and theoretical predictions, in 
form of the wind-momentum luminosity diagram (left panel: 
"Method 1"; right panel: "Method 2"). Circles correspond to 
our investigation, denoted as open for objects with f c \ = 1, 
and filled for those two objects which show clear evidence for 
clumping. 

The positions found for our objects scatter in a similar way 
around the "observed" WLR for Galact ic O-type supergiants 



and giants/dwarfs (dotted/dashed, from iMarkova et al. 2004) 



and the corre sponding theoretical prediction s (solid, from 



Vink et al. 2000l) as the results obtained by Kudritzki et al 



(2006, asterisks) for a different sample, with He 2 - 108 in com- 
mon. We also show the results from IPauldrach et"al] (2004, 
crosses), obtained from a UV-analysis based on a consistent 
solution of the hydrodynamics. 

Leaving aside the most extreme outliefl those two CSPN 
with clumping corrections lie close to but below the O-star 
WLR, one object is on the WLR and the other two well above. 
Since we are not able to exclude clumping corrections for the 
latter objects, it might b e that they also follow the O-star trend. 



As already outlined by lKudritzki et al.l (|2006), however, a def- 



inite statement regarding the question of whether CSPN are a 
lo w luminosity extension of the O-star relation (a s suggested 
by IPauldrach et al.|[l988l and iKudritzki etaill 19971) or whether 
there is no such clear relation (see also Tinkler & Lamersl2002 ) 
is still not possible at the present stage. 



8 corresponding to M 1-37 (uppermost left circle), which is also 
the most extreme object in our sample (with very strong Sim lines) 
and might be in the transition phase to a WR type object, in which 
case it must have a larger M and wind-momentum rate than "normal" 
CSPN. 
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Fig. 13. Modified wind-momentum rates, D mom , as a function of luminosity. Left: results from "Method 1"; right: results from 
"Method 2". Circles correspond to the results for our sample (open: f c \ = 1, filled: f c \ + 1), asterisks to the results from 
Kudritzki et al. (2006) and pluses to the results from lPauldrach et all ( 2004 ). The s pecial symbols den ote the various solutions for 
He 2-108, as in Figure[TT] The solid line displays the theoretical prediction from Vink et al. ( 2000l). and the dotted and dashed 
line the "observed" WLR for Galactic O-type supergiants and giants/dwarfs, respectively (from lMarkova et aTll2004l) . The object 
that deviates most from the predicted WLR is the most extreme object in our sample, M 1-37. 



There is one additional argument which might support the 
former hypothesis. If one claims that the CSPN indeed form 
such an extension, one obtains a certain constraint on their 
wind-strength parameter, log Q: Rewritten in terms of this pa- 
rameter, the WLR for a certain class of objects reads 

log£> mom » x\ogL + D, i.e., 

log Q + 2 log/?. + - log Voo * 4 x log r eff + 2 x log R m + D' 



(Voo <* v esc = y/gRt) 
loge + (2-2x+^)log/?, « 4xlog7; ff -^logg + D", (5) 

with x and D the slope and offset of the specific WLR. If now 
the CSPN would follow such a WLR and were an extension of 
the O-star relation, this equation implies 



(log Q - 0.5 log/?»)cspN = (log Q - 0.5 log/?*) 

O- stars » 



(6) 



if we compare CSPN and O-stars of similar effective temper- 
atures and gravities, using a typical slope of l/x = 0.53 for 
the O-star WLR. In other words, the average wind-strengths of 
CSPN must be roughly 0.5 dex larger than the corresponding 
O-star values at same T e ff and log g, if we adopt a typical dif- 
ference of 1 dex in radii. But this is what we actually see in our 
spectra (cf. Section [3~5l i: Compared to O-supergiants of similar 
spectral type, all our sample stars show much more emission 
in Ha, immediately pointing to a higher log Q. Indeed, the av- 
erage value for our sample is log Q = -11.9, compared with 
typical (late) O-supergiant values in between -12.8 to -12.4. 

Summarizing, a similar WLR implies that for CSPN one 
expects larger Q (for smaller /?» and v m ), whereas, vice versa, 
the observed, strong Ha emission might imply that the CSPN 
indeed lie on the O-star WLR. Alternatively, there is the possi- 
bility that the clumping properties of O-stars and CSPN differ 
significantly, with stronger clumping in CSPN. 



5. Discussion 

This project has produced a remarkable mixture of encourag- 
ing and alarming results. On the positive side, leaving aside 
M 1-37 with its unique spectrum and uncertain \ogg, the four 
remaining bulge central stars turn out to have "Method 1" spec- 
troscopic distances within +10% to +15% (depending on the 
extinction law) of their average distance, well within the ex- 
pected uncertainties. The alarming fact is of course that their 
average d istance, if we adopt t he reddening corrections sug- 
gested by lStasinska et alJ (119921) . is a factor 1.5 larger than the 
currently accepted distance to the Galactic bulge. This allows 
the following possible interpretations: 

(1) If the reddening corrections based on a low /?y and radio- 
H/3 extinctions are correct, then the combination of our 
model atmospheres plus standard post-AGB evolutionary 
tracks gives systematically too large spectroscopic dis- 
tances, despite their excellent internal agreement. 

(2) If the standard reddening law is valid, and using C opt , then 
our "Method 1" spectroscopic distances give the right dis- 
tance to the Galactic bulge to within 13%, which would 
have to be considered a successful test; but then we need 
to find a reason for the observed systematic difference be- 
tween the optical and radio-H/3 extinctions. 

In the present situation we do not feel that we have enough 
information to decide between (1) and (2). Our referee, A. 
Zijlstra, has pointed out that our sample objects have (by bad 
luck) relatively poorly known radio fluxes. Most come from 
very old radio data, obtained at the very start of the VLA, and 
we cannot be as confident of these as we are of more recent ob- 
servations. Only He 2-260 can be considered as reliable; the 
others need confirmation. 

The analysis of our comparison object from the "solar 
neighborhood", He 2-108, reveals a prototypical CSPN mass, 
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of about 0.6M Q (both for the low and high T e ff solution) if 

"Method 1" is followed. This is in agreement with the 0.67M o 

I 1 J L 

derived by Mendez et al. ( 1992) and the recent value of 0.63M G 

by iKudritzki et alj (12006I) . The higher masses derived using 
"Method 2", however, lacks any significance (except for the 
fact that it overlaps within the error bars with our "Method 1"- 
mass), due to its highly uncertain distance (adopted from simi- 
lar analys es) and the large er ror propagation of this quantity. 
Recently, Napiwotzki d2006l) has claimed, using an analysis 
of kinematical parameters based on simulated Galactic orbits, 
i.e., by means of a completely different approach, that a low 
mass similar to our "Method 1" result would be more compat- 
ible with the kinematical parameters one would expect from a 

member of the thin-disk population. 

r — — 1 J k 

On the other hand, Pauldrachet al. (2004) derived a sig- 
nificantly larger mass of 1.33M for this object. In contrast to 
our approach, they have solved the hydrodynamical equations 
in a self-consistent way, and compared the corresponding, syn- 
thetic UV spectra with observations. T e g then follows from the 
observed ionization equilibrium, whereas mass and radius are 
uniquely constrained from the terminal velocity (being propor- 
tional to the photospheric escape velocity, oc (M,//?,) 1 ^ 2 , and 
mass-loss rate (depending on luminosity, i.e., again on stellar 
radius). Note that both wind parameters are clearly "visible" in 
the UV-lines. 

Ideally, the results of optical and UV spectroscopy should 
be similar. By comparing the derived effective temperatures 
and gravities for our high T e ff solution (Fig.QTl 'H' vs. 'P'), at 
least for these parameters this is indeed the case (i.e., the ratios 
M*/Rl do agree). Thus, the difference in mass is due to a differ- 
e nt radius (implying diff erent ratios M*/R*) which in the case 
of lPauldrach et alj d2004l) results from hydrodynamical consid- 
erations, whereas here we have derived the mass directly from 
evolutionary tracks. 

A solution of this discrepancy cannot yet been provided, 
but will be commented on in the following. At first, let us 
give a few additional remarks con cerning PN distances. I n the 
same publication as just referred to. lPauldrach et alJ (|2004) pre- 
sented a comparison between spectroscopic and other PN dis- 
tance determinations, which will not be repeated here except 
to comment on the fact that one of the very few cases of a 
clear discrepancy (Hipparcos parallax of PHL 93 2) has been 
resolv ed in favor of the spectroscopic distance. Harris et al] 
(2007) reports new trigonometric parallaxes for several cen- 
tral stars, pointing out that the Hipparcos parallax of PHL 932 
is wrong. The new geometric distance of PHL 932 is 300 pc, 
in perfect agreement with the spectroscopic distance we ob- 
tain if we assume a stellar mass of 0.3 M and the follow- 
ing atmospheric parameters: T e g = 36 000 K and logg = 5.7, 
which a re the averages from two independent spec troscopic 
studies (Men dez et alJll988al and lNapiwotzkil [l999). Clearly 
the Harris geometric distance confirms the small stellar mass 
and proves definitely that PHL 932 cannot be a post-AGB 



star (Mendez et al. 1988a). So here we have a case where 



Of course we must add that those central stars with trigono- 
metric parallaxes are high-gravity stars with no wind features; 
so this agreement is not a confirmation of the wind models. The 
same comment applies to the central star of K 648, the PN in 
the globular cluster M 15, whose spectroscopic distance is in 



excel lent agreement with the cluster distance (IMcCarthy et al 



1997). Let us argue in the following way: if the combination of 



the model atmospheres pass the test. lHarris et al ] (E007n also 
shows agreement, within the uncertainties, between his paral- 
laxes and spectros copic distances for several PN, mostly from 
iNapiwotzkild 1999b . 



post-AGB theory and windless model atmospheres gives good 
distances (as indicated by the high-gravity stars and K 648) 
and the combination of post-AGB theory and wind models 
fails systematically, then it would be easy to conclude that 
the wind models must be wrong. But this is not the case - 
evidence for such systematical failures have not been found 
so far, and it should be noted that our Galactic bulge sam- 
ple includes central stars with weak and strong winds, yet 
all their "Method 1" spectroscopic distances agree, with no 
differences that could be attributed to wind strength effects. 
Alternatively, we have the possibility that another group of 
CSPN with w inds exists (e.g., so me objects from the sample 
analyzed by iPauldrach et alj 120041) . which follow a different 
evolutionary path and for which the standard CSPN tracks do 
not apply. A further indi cation of such a poss ible qualitative 
difference is the fact that Kudritzki et al. (2006), using a simi- 
lar method as in this paper, also based on the standard evolu- 
tionary tracks, found higher masses ( around 0.74 M ) for th e 
same sample of objects as analyzed by IPauldrach et al.1 d2004l) . 
whereas the mean mass of the CSPN in our current sample is 
around 0.68 M Q . 

In summary, the combination of excellent internal agree- 
ment (and some additional supporting evidence) plus the 
extinction problem leads us to conclude that our Galactic bulge 
distance test is undecided. 

We can see two ways of trying to force a decision: 

(1) a careful, multi- wavelength redetermination of the inter- 
stellar extinction toward the bulge PN we have studied, 

(2) enlarge the sample of bulge central stars, trying to get high- 
quality, high-resolution spectrograms with minimum nebu- 
lar contamination. 

Alternatively, as already emphasized in the introduction, post- 
AGB tracks and model atmospheres could be independently 
tested as soon as medium-resolution spectroscopy of CSPN in 
the Magellanic Clouds, with efficient nebular light subtraction, 
becomes feasible; the advantage being the small amount of red- 
dening in that direction. 

6. Summary and conclusions 

We have selected a small sample of PN central stars in the 
Galactic bulge and obtained high-resolution spectrograms with 
the Keck HIRES spectrograph. We present spectral classifica- 
tion of these central stars and we model their optical spectral 
features (including wind effects) using a state of the art non- 
LTE, spherically symmetric model atmosphere code that in- 
cludes the effects of line blocking/blanketing and wind clump- 
ing. We obtain effective temperatures, surface gravities, wind 
parameters and He and Si abundances. With one exception, He 
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and Si abundances are normal for H-rich central stars. The ex- 
ception is M 1 -37, a central star with spectacular P Cygni pro- 
files indicative of a very dense wind, He-rich, and with an ap- 
parently very high Si abundance, 28 times solar, unexplained 
for the moment. This star also shows strong unidentified emis- 
sions at 4485 and 4504 A, which we have refrained from iden- 
tifying with Silicon. 

Armed with the stellar atmospheric parameters we then cal- 
culate other parameters not directly derivable from our analy- 
sis, using two methods. "Method 1" uses post-AGB evolution- 
ary tracks to read masses, from which stellar radii, luminosities, 
mass loss rates and spectroscopic distances follow. "Method 2" 
assumes the distance to the Galactic bulge and derives stellar 
radii and other quantities; this method turns out to be very un- 
certain because of the strong dependence of mass on stellar ra- 
dius, which is itself not well determined. 

If we focus our attention on the results of "Method 1", the 
most remarkable conclusion regards the internal agreement of 
the spectroscopic distances of at least four bulge central stars 
within our sample (leaving aside the extreme object M 1-37 
with its ill-defined gravity): all of them are quite the same to 
within +10% to +15% (depending on the extinction law) of 
their average distance. This result is well within the expected 
uncertainties. There is a problem, however: the average dis- 
tance to the Galactic bulge is strongly dependent on the adopted 
reddening correction, to such an extent that we cannot decide 
if "Method 1" passes the consistency test for our sample (if we 
adopt the most recently derived reddening corrections, the aver- 
age distance turns out to be a factor 1 .5 larger than the currently 
accepted distance to the Galactic bulge). Additional studies of 
extinction for the PN in our sample will be necessary, and it 
would be good to increase the size of the bulge CSPN sample 
subject to this kind of analysis, while we wait for a chance to 
repeat this kind of study in the Magellanic Clouds, which are 
almost free of the extinction problem. 
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